Signatures of the superfluid-insulator phase transition in laser driven dissipative 

nonhnear cavity arrays 



A. Tomadin/ V. Giovannetti/ R. Fazio, ^ D. Gerace,^ I. Carusotto,^'^ H.E. Tiireci,^ and A. Imamoglu^ 

^NEST-CNR-INFM and Scuola Normale Superiore, Piazza dei Cavalieri 7, 1-56126 Pisa, Italy 
^CNISM and Dipartimento di Fisica "A. Volta," Universitd di Pavia, 27100 Pavia, Italy 

^Institute for Quantum Electronics, ETH Zurich, 8093 Zurich, Switzerland 
^BEC-CNR-INFM and Dipartimento di Fisica, Universitd di Trento, 1-38050 Povo, Italy 

(Dated: June 15, 2010) 

We analyze the non-equilibrium dynamics of a gas of interacting photons in an array of coupled 
dissipative nonlinear cavities when driven by a pulsed external coherent field. Using a mean- field 
approach, we show that the response of the system is strongly sensitive to the underlying (equilib- 
rium) quantum phase transition from a Mott insulator to a superfluid state at commensurate filling. 
We find that the coherence of the cavity emission after a quantum quench can be used to determine 
the phase diagram of an optical many-body system even in the presence of dissipation. 
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Recent theoretical advances in cavity Quantum Elec- 
trodynamics (QED) have indicated arrays of coupled 
nonlinear cavities as potential candidates to explore 
quantum many-body phenomena of light [1 . Initial pro- 
posals for realizing a Mott phase of polaritons [2^ have 
been scrutinized in great detail [3 and different schemes 
have been put forward to simulate a variety of correlated 
quantum states [4^7]. Given the importance of dissipa- 
tion in state-of-the-art solid-state (single) cavity QED 
devices (see for example [8 ), experiments in cavity QED 
arrays are expected to be performed under strongly non- 
equilibrium conditions, with an external source compen- 
sating for the loss of photons (see Fig. [T]). These condi- 
tions, dictated by the experimental constraints, put for- 
ward QED arrays as natural candidates to explore the 
rich world of non-equilibrium quantum many-body sys- 
tems [10 -.12 , that is the subject of recent interest in the 
context of ultracold atoms as well [I3l [14] . At the same 
time new important questions arise related to whether 
or not it is possible to realize and detect, under realistic 
non-equilibrium conditions, the very rich phase diagram 
[1] predicted for QED arrays at equilibrium. 

In this Letter we propose an answer to this last ques- 
tion thus providing the missing link between the initial 
ideas of cavity arrays as quantum simulators and future 
experimental realizations. Specifically we will show how 
to detect the superfluid-insulator phase transition for a 
dissipative nonlinear cavity array driven by an external 
pulsed laser. In analogy to what happens in a quan- 
tum quench [151 HB] ^ photons are excited in the cavities 
by a periodic train of short, coherent pulses, and then 
evolve according to the complex many-body dynamics de- 
termined by the simultaneous action of hopping, strong 
interactions and losses. The properties of the secondary 
emission are measured in the time lapse between subse- 
quent pulses: even though the losses drive the system to- 
wards the vacuum state, distinct signatures of the quan- 
tum phase transition are found in the statistics of the 




FIG. 1: (Color online) Schematic representation of a possible 
realization of the proposed system, employing a 2-dimensional 
array of photonic crystal defect cavities [9 . Photons in the 
cavities have a finite lifetime, which requires losses to be 
compensated by some external pump, in our case a coher- 
ent pulsed laser field. Nonlinear processes in each cavity lead 
to an effective Bose-Hubbard repulsion between the cavity 
photons 1 . Neighboring cavities are tunnel coupled to each 
other by spatial overlap of the cavity mode profiles. 

transient emitted light. First- and second-order photon 
correlation are shown to be powerful tools to detect the 
existence of a non-zero order parameter as they are very 
sensitive to the delocalization of photons through the ar- 
ray, although more complete measurement schemes exist 
[17J. 

The time evolution of the array in the presence of driv- 
ing and losses is described by the master equation 

dtp{t) = -i[n,p] + c[p] , (1) 

for the density matrix p of the system. The first term 
on the r.h.s. of Eq. ([T]) accounts for the unitary evolu- 
tion of the system while the second term accounts for 
the damping. As discussed in [1 , an array of coupled 



cavities in the absence of losses can realize an effective 
Bose-Hubbard model whose Hamiltonian reads 

7i = - J ^ aja^/ + ^ X] - 1) 
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+ A^n, + ^[s(t)a, + £(t)*aJ] , (2) 
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with a^, aj, and hi denoting, respectively, the annihila- 
tion, creation and number operator of the electromag- 
netic mode of the ith cavity. The first term of Eq. ([2|, 
quantified by the coupling constant J, describes the hop- 
ping of photons between neighboring cavities (the lat- 
tice has a coordination number z). The second term, 
quantified by [/, represents the effective nonlinearity aris- 
ing from the specific light-matter interaction mechanism. 
The achievable values of J and U depend on the specific 
implementation of the model [1 . The last two terms in 
Eq. ([2| describe the coupling of the photons in the cavi- 
ties with a uniform coherent pump in the rotating frame 
- here €{t) is the (slowly- varying) envelope of the external 
driving field and A is the detuning of the laser frequency 
from the bare cavity resonance. The dominant source of 
dissipation in the system is the leakage of photons out 
of the cavities. In the Markov approximation, the de- 
cay process is described by a Liouvillian in the Lindblad 
form C[p] = i^'^ii'^CLipal — hip — phi)^ where is 
the photon lifetime. 

Solving Eq. ([T]) is a formidable task, as it describes 
an open quantum many-body system out of equilibrium. 
We thus resort to a self-consistent cluster mean-field ap- 
proach, which should give reliable results if the array is 
(at least) two-dimensional. The cluster mean-field (for 
its equilibrium version see for example [18 ) is based on 
approximating the Hamiltonian ([2| with that of a clus- 
ter of cavities that interact with the rest of the lattice 
via a mean-field term, i.e. Hmf = reQc^^^^' ~ 
'^^{t/yieni./'^nc^^l^^'^'^'^']' Here, is the Hamil- 
tonian PI restricted to the sites within the cluster (ftc)- 
The second term instead is the mean-field expression for 
the hopping from the cavities £ in the cluster boundary 
(l^b) to their nearest-neighbors £^ outside the cluster. It 
is a function of the time-dependent mean-field order pa- 
rameter tlJi{t) which, for the ^th cavity, is determined 
by the self-consistency condition ^^^(t) = Ti [a i p{t)]. In 
the remainder of the paper we will consider only the two 
cases in which the cluster consists of one or two sites. 
For the numerical integration of the master equation we 
use a fourth order Runge-Kutta algorithm and truncate 
the local Fock basis {\n)i}]^'^Q' for the ^th cavity, with an 
upper cutoff n^ax — 20. 

In order to simplify the presentation, we first consider 
an idealized case in which the laser pulse sets, at t = 0, 
the cavity array in a Fock state pin = with 
one photon per cavity. We will discuss the precise shape 
of the required pulses and the robustness of the response 
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FIG. 2: Time evolution of the zero-time delay second-order 
correlation function and of the rescaled order parameter (in- 
set), after a pulse that projects the cavities into a Fock state 
with one photon per cavity. The rescaled order parameter 
'0(t) = \^lj{t)\/ ^/n(t) compensates for the decay of the photon 
number in the cavities and is related to the coherent fraction 
of the emission {k = 10"^ U, zJ = 3.0 U, A=0). 

of the system to imperfections in the pulse shape further 
below. We first focus on single-site mean field and an- 
alyze the coherence properties of the light emitted from 
the cavity array after each pulse. The time-dependence 
of the photon population n{t) of each cavity is trivial, 
the corresponding equation of motion can be integrated 
yielding an exponential decay n{t) = n(0)e~^'^^. Much 
more interesting are the properties of the order parame- 
ter '0^(t) and the zero-time delay second-order correlation 
function, g2 = {hj — h^) / {h^)'^ ^ for t > 0, i.e. after the 
laser pulse is switched off. If J /U < {J/U)c {{J/U)c is 
the value at which the transition from the superfluid to 
the Mott state occurs) any initial fluctuation in the order 
parameter does not grow in time and the density matrix 
in photon-number representation is essentially diagonal 
throughout the relaxation process. In the opposite case, 
J/U > (J//7)c, a superfluid instability can develop and 
leads to nonlinear oscillations of both |'0(t)| and ^2- This 
kind of instability has been discussed, in the absence of 
dissipation, in [15 . Here we show that the instability is 
present also in the presence of losses. It can be gener- 
ated by means of pulsed laser and detected by measuring 
the properties of the light emitted from the cavities. On 
the contrary, within our mean- field analysis, continuous 
shining of coherent light on the cavities would wash out 
the effects of the instability. 

For the open system we are considering, a generic ex- 
ample of the instability induced by the pulsed pump is 
shown in Fig. |2j where we show the evolution of g2{t) 
and (normalized) \^p{t) \ calculated for t > 0. The dynam- 
ics is characterized by three different regimes. At short 
times, a linear instability sets in and both quantities in- 
crease exponentially, notwithstanding the slow decay of 
the photon population. At intermediate times, the collec- 
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tive nonlinear dynamics of the array leads to oscillations. 
In the long time limit, these oscillations damp out with 
a time constant of the order of k.~^ (in this regime, when 
n{t) <C 1, it can be shown that \ip{t)\ ~ ^^nit)). 

An experimentally measurable observable that yields 
a clear signature of the different regimes in the initial 
transient is the zero-time delay second-order correlation 
function averaged over a certain interval of time. To cal- 
culate such time averaged observables, we integrate the 
equation of motion Eq. ^ up to the time 2.0 a^"^, and de- 
termine the time averaged values {\ip{t)\)t and {g2)t in the 
time interval 1.0 a^"^ < t < 2.0 a^"^. Considering time- 
integrated quantities allows to relax the experimental re- 
quirements on the time resolution of the measurements, 
which can hardly exceed in realistic state-of-the-art 
devices. 

In Fig. [3] we show (|?/^(t)|)t and {g2)t as a function of 
the ratio J/U^ for different values of the dissipation rate. 
Both quantities vanish if J/U is smaller than the critical 
value, meaning that the order parameter does not develop 
any instability, and consequently the light emitted from 
the cavities is strongly antibunched. In the opposite case 
both {\tlj{t)\)t and {g2)t are different from zero. The time- 
averaged g2 as a function of J/U monotonically rises from 
zero to almost unity, thus showing a crossover from an- 
tibunched to Poissonian statistics. Integrating the equa- 
tion of motion for times larger than 2.0 leaves the 
features depicted in Fig. |3] qualitatively unaltered. We 
emphasize that the crossover of both {\ilj{t)\)t and {g2)t 
takes place in a very narrow range of J/U values around 
the transition point of the closed system. The sensitivity 
of the light statistics to the coupling between neighbor- 
ing cavities is a signature of the many-body origin of this 
phenomenon. The second-order correlation is thus an ex- 
cellent candidate to detect the different quantum phases 
in cavity arrays. 

A linear stability analysis, valid for t [/, t <C 1, shows 
that the rate of the initial exponential build-up of the or- 
der parameter is reduced by a quantity hc^ at fixed J/U] 
for large enough the superfluid instability is then en- 
tirely suppressed. As a consequence, the critical ratio 
{J/U)c increases by /{2U'^). Such trend is already 
visible for U/n 10^. We note that this suppression of 
the superfluid phase is not a consequence of the coupling 
to a bath [13]. It is important to note that U / n c::^ 10^ is 
experimentally achievable in a solid-state cavity with a 
linewidth of :^ 10 /ieV (corresponding to Q 10^ at opti- 
cal or near-infrared frequencies) and interaction strength 
~ 1 meV, which are within reach in current solid-state 
cavity QED systems [1£ . In the case of coupled pho- 
tonic crystal cavities (see Fig. [T]) the range of param- 
eters that one can expect (assuming = 0.1 meV) is 
< J/U < 50 [20]. Alternatively, circuit-QED also pro- 
vides an ideal experimental realization for the present 
proposal [2T. 

All the results presented here are robust even when 
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FIG. 3: Time- average of g2{t) (main panel) and integral of 
|'0(t)| (inset), in the time interval 1.0 < t < 2.0 Hi~^ for 
different values of the damping: k, = 2.0 x 10~^ U (empty 
circles), 1.0 x 10"^ U (filled circles), and 0.5 x 10"^ U (empty 
triangles); A = 0. The vertical dashed line marks the critical 
value of the Mott-to-superfluid transition for the equilibrium 
Bose-Hubbard model {{zJ/U)c ^ 0.17). 



some of the assumptions made so far are relaxed. In 
the vicinity of the critical point, the order parameter 
evolves on time scales which are much longer than the 
cavity parameters [15 . In contrast, the correlation func- 
tion might show a more complex short-distance evolution 
which could invalidate the analysis presented above. In 
order to check this possibility, we carried out the same 
analysis in the two-site mean field approximation, where 
both the dynamics of the order parameter and the short 
time scale dynamics governed by J and U are present. All 
the results shown so far are fully confirmed. In addition, 
we checked that the uniform solution for the order param- 
eter is stable against possible spatial inhomogeneities of 
the initial state. 

We also checked that our findings do not depend on the 
exact form of the state after the pulsed excitation. We 
considered an initial state that is not necessarily pure, 
but that entirely projects the system onto the subspace 
spanned only by the two Fock states \0)£ and |1)^. We 
considered depletion of the average filling up to 20% and 
scanned the whole range of initial coherences. For all the 
possible initial states, {g2)t in the superfluid regime was 
always markedly different (at least three orders of mag- 
nitude) from the insulating one, as shown in Fig. |4] This 
insensitivity to the initial conditions in the superfluid 
regime is due to the nonlinearity of the time-evolution 
for times shorter than the photon lifetime. In fact, in the 
insulating regime the initial correlations are suppressed 
due to the photon blockade, while in the superfluid re- 
gion even the absence of any initial correlation is quickly 
compensated by the cooperative action of the photons. 
We remark that while {g2)t predicts strongly antibunched 
cavity emission in the insulating regime for all considered 
variations in the initial states, the actual degree of anti- 
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FIG. 4: Time average of g2{t), for n = 10~^U. The pop- 
ulation po of the vacuum in the initial state is po = 0.02 
for the filled symbols and po = 0.2 for the empty symbols. 
{zJ/U = 0.1 (triangles) and zJ/U =1.0 (circles), below and 
above the critical value; A = 0.) 



FIG. 5: Optimized pulse envelopes for the preparation of a 
Fock state with one photon per cavity, with a fidelity T ^ 
99%. The dashed line is the pulse e^^\t) for J = 0. The 
solid line is the pulse e^'^\t), detuned by A = 43.25 At, with 
zJ/U = 0.5. The duration Tp = 0.1/ k of the pulse is much 
shorter than the photon lifetime — 10~^ U in the plot). 



bunching does depend on the initial state coherence; at 
present, we do not understand this feature. 

We now discuss pulse shapes that project the coupled 
cavity array into a desired initial state (such as pin dis- 
cussed before). A 7r-pulse applied to an isolated cavity, 
drives the vacuum state into the desired initial state, pro- 
vided that the duration of the pulse is shorter than the 
photon lifetime i<i~^ . To design a suitable pulse shape 
for the whole array, we resorted to a rapidly-converging 
quantum optimal control algorithm [2^: first, we found 
an optimized envelope e^^\t) (see Fig.^ dashed line) for 
a single cavity (J = 0). As dissipation does not play a 
substantial role during the application of the pulse in the 
time interval — Tp < t < 0, we use an algorithm of quan- 
tum optimal control for pure states, obtaining a fidelity 
99% with respect of the desired state pin- In the 
case of an array of coupled cavities, J > substantially 
modifies the response of the system to the pumping field 
and it is necessary to devise a different shape e'^^^t) of 
the pulse for each value of the hopping amplitude (see 
Fig. [sj solid line). It is nevertheless possible to prepare 
the initial state with high fidelity (99% in Fig. [5| in a 
substantial range of values up to zJ ^ U. 

In conclusion, we have shown that the coherence prop- 
erties of the secondary cavity emission, induced by a 
pulsed excitation, provide signatures of the collective 
many-body phase of the array of nonlinear cavities. Our 
analysis fully takes into account the intrinsically driven- 
dissipative nature of the strongly correlated quantum 
many-body system and identifies how dissipation influ- 
ences the underlying equilibrium quantum phase transi- 
tion that this system exhibits in the absence of losses. 
Coupled cavities are promising candidates to simulate 
open quantum many-body systems and a study of their 
non-equilibrium behavior will certainly unveil a variety 
of new phenomena. 



We acknowledge useful discussions with Sebastian 
Schmidt. This work was supported by EU IP-EuroSQIP 
(RF), the Swiss NSF under Grant No. PP00P2-123519/1 
(HET) and an ERG Advanced Investigator Grant (AI). 



[1] M.J. Hartmann, F.G.S.L. Brandao, and M.B. Plenio, 
Laser & Photon. Rev. 2, 527 (2008). 

[2] M.J. Hartmann, F.G.S.L. Brandao, and M.B. Plenio, Na- 
ture Phys. 2, 849 (2006); A.D. Greentree et al, Nature 
Phys. 2, 856 (2006); D.G. Angelakis, M.F. Santos, and 
S. Bose, Phys. Rev. A 76, 031805(R) (2007). 

[3] D. Rossini and R. Fazio, Phys. Rev. Lett. 99, 186401 
(2007); N. Na et a/., Phys. Rev. A 77, 031803(R) (2008); 
M. Aichhorn et aL, Phys. Rev. Lett. 100, 216401 (2008); 
S. Schmidt and G. Blatter, ibid. 103, 086403 (2009). 

[4] M.J. Hartmann, F.G.S.L. Brandao, and M.B. Plenio, 
Phys. Rev. Lett. 99, 160501 (2007). 

[5] A. Kay and D.G. Angelakis, EPL. 84, 20001 (2008). 

[6] J. Gho, D.G. Angelakis, and S. Bose, Phys. Rev. Lett. 
101, 246809 (2008). 

[7] I. Garusotto et a/., Phys. Rev. Lett. 103, 033601 (2009). 

[8] K. Hennessy et a/.. Nature 445, 896 (2007); A. Faraon 
et al, Nature Phys. 4, 859 (2008). 

[9] H. Altug and J. Vuckovic, Appl. Phys. Lett. 84, 161 
(2004). 

[10] D.E. Ghang et a/.. Nature Phys. 4, 884 (2008). 

[11] D. Sarchi et a/., Phys. Rev. B 77, 125324 (2008); D. Ger- 
ace, et al, Nature Phys. 5, 281 (2009). 

[12] I. Garusotto and G. Giuti, Phys. Rev. Lett. 93, 166401 
(2004); M.H. Szymahska, J. Keeling, P. B. Littlewood, 
ibid. 96, 230602 (2006); M. Wouters and I. Garusotto, 
ibid. 99, 140402 (2007); A. Amo et a/., Nat. Phys. (in 
press) . 

[13] S. Diehl et a/.. Nature Phys. 4, 878 (2008). 

[14] A.J. Daley et al, Phys. Rev. Lett. 102, 040402 (2009). 

[15] E. Altman and A. Auerbach, Phys. Rev. Lett. 89, 250404 



5 



(2002). 

[16] R Calabrese and J. Cardy, Phys. Rev. Lett. 96, 136801 
(2006); C. Kollath, A.M. Lauchli, and E. Altman, ibid. 
98, 180601 (2007); M. Rigol, V. Dunjko, and M. Olshanii, 
Nature 452, 854 (2008). 

[17] M.F. Santos et al, Phys. Rev. A 63, 033813 (2001). 

[18] T. Oguchi, Prog. Theor. Phys. 13, 148 (1955). 

[19] S. Combrie et a/., Opt. Lett. 33, 1908 (2008); A. Verger 



et a/., Phys. Rev. B 73, 193306 (2006). 
[20] S. Vignohni et a/., Appl. Phys. Lett. 94, 151103 (2009); 

S. Vignohni et al, ibid. 96, 141114 (2010). 
[21] D.L Schuster et al, Nature 445, 515 (2007); L.S. Bishop 

et a/., Nature Phys. 5, 105 (2009). 
[22] J. Werschnik and E.K.U. Gross, J. Phys. B 40, R175 

(2007). 



